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Introduction 

Nature  of  the  Problem 

While  it  is  hoped  that  molecular  detection  and  intervention  will  one  day  provide  a 
more  effective  treatment  modality,  currently,  in  the  words  of  the  Report  to  the  U.S.  Army 
Medical  Research  and  Development  Command  on  the  Strategies  for  Managing  the 
Breast  Cancer  Research  Program,  "no  dominant  etiology  for  breast  cancer  has 
emerged. ..[this]  would  lessen  the  prospects  for  any  quick  and  easy  prevention 
strategies....  mammography  is  the  method  of  choice  for  screening  women  to  detect 
breast  cancer ...  mammography  has  proven  to  be  the  most  effective  means  of  reducing 
breast  cancer  morbidity  and  mortality." 

The  primary  theoretical  limitations  of  mammography  are  the  system  resolution, 
which  determines  the  minimum  size  of  the  detectable  malignancy,  and  the  need  to 
expose  the  patient  to  ionizing  radiation.  Using  innovative  new  technology  to  improve 
system  resolution  and  reduce  required  dose  will  increase  the  effectiveness  of  this 
proven  screening  modality,  with  a  direct  and  immediate  impact  on  mortality.  In  addition, 
this  direct  digital  system  can  avail  itself  of  the  advantages  of  digital  processing, 
including  improved  image  contrast  and  resolution  at  reduced  radiation  dose.1  In 
practice,  mammographic  imaging  is  often  limited  by  quality  assurance  issues,  which 
can  also  be  favorably  addressed  by  digital  processing. 

Purpose 

The  overall  objective  of  this  proposal  is  to  develop  a  mammographic  system  with 
extremely  high  scatter  rejection  and  dynamic  range,  good  resolution  and  low  patient 
dose.  This  will  be  accomplished  by  developing  a  direct  x-ray  detector  interfaced  with  a 
capillary  x-ray  optic  in  an  appropriately  designed  mammographic  system. 

Kumakhov  capillary  x-ray  optics,  invented  in  the  mid-eighties,  provide  an 
innovative  new  way  to  control  x-ray  beams.  Such  optics  will  provide  extremely  efficient 
scatter  rejection,  while  allowing  beam  magnification,  demagnification,  and  shaping  to 
match  with  the  newly  developing  high  efficiency  direct  x-ray  detectors.  These  detectors 
owe  their  high  efficiency  and  resolution  to  the  direct  detection  of  x-ray  photons  without 
requiring  phosphors  for  the  conversion  to  visible  light.  An  integrated  system  of  optics 
and  detectors  will  be  developed  in  a  highly  collaborative  effort  involving  recognized 
leaders  in  the  fields  of  capillary  optics,  x-ray  detectors,  digital  radiology,  and 
mammography.  Testing  will  be  performed  on  each  of  the  elements  independently,  and 
as  an  integrated  unit  in  a  mammographic  system. 
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Background 

Benefits  of  Digital  Mammography 

Conventional  film/screen  mammography  suffers 
from  limited  dynamic  range  and  film  granularity,  which  can 
reduce  the  sensitivity  of  detection  of  microcalcifications. 

Digital  detection  provides  high  dynamic  range,  improving 
contrast,  and  greatly  increasing  the  tolerance  of  the  final 
image  to  under-  or  over-exposure.  Digital  images  can  be 
enhanced  and  are  amenable  to  computer  aided  diagnosis. 

Spectral  information  can  be  included  if  it  is  available. 

Finally,  digital  images  can  be  quickly  transported  for  skilled 
consultation. 

Capillary  Optics 

Kumakhov  capillary  optics  are  bundles  of  hollow 
glass  capillary  tubes  with  inner  diameters  as  small  as  a  few 
microns.  A  typical  fiber  is  shown  in  Figure  1 .  X  rays 
incident  on  the  interior  of  the  glass  tubes  at  small  angles 
can  be  guided  down  the  tubes  by  total  external  reflection. 

The  polycapillary  fibers  guide  x  rays  in  a  manner  analogous  to  the  way  fiber  optics 
guide  light.  Arrays  of  curved  tapered  fibers  can  be  used  to  focus,  collimate  and  filter 
x-ray  radiation.2’3-4’5-6  Such  arrays  can  be  manufactured  by  stringing  hollow  glass 
polycapillaries  through  metal  grids,  as  shown  in  Figure  24,  or  manufactured  without 
grids  as  a  monolithic  optic,  as  shown  in  Figure  2. 

The  critical  angle  for  total  external  reflection  of  x  rays  by  glass  polycapillaries  is 


where  cop  is  the  plasma  frequency  of  the  glass,  about  30  eV,  and  to  is  the  photon 

frequency.  The  critical  angle  is  1.5  milliradians  at  20  keV.  The  x  rays  can  be 
transmitted  in  a  curved  tube  so  long  as  the  tube  is  small  enough  and  bent  gently 
enough  that  the  angles  of  incidence  are  kept  less  than  the  critical  angle.  For  a  given 
radius  of  curvature,  this  requires  increasingly  small  diameter  tubes  as  the  x-ray  energy 
is  raised.  In  order  to  avoid  the  mechanical  limitations  of  such  small  tube  sizes, 
polycapillary  fibers  are  employed  with  channel  sizes  (typically  1-30  ^m)  much  smaller 

than  the  outer  diameter  (300-1000  p.m). 

The  use  of  Kumakhov  capillary  optics 
in  place  of  more  conventional  scatter 
reduction  grids  in  a  mammographic  system 
has  significant  potential  to  provide  improved 
resolution,  increased  contrast 
enhancement,  and  reduced  dose  in 
mammographic  imaging.  The  optics  can 
also  be  used  to  mate  the  radiographic 
image  with  a  digital  detector  by  appropriate 


Figure  2 .  Monolithic  optic  geometry. 
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Figure  1.  Cross  sectional 
scanning  electron  micrograph 
of polycapillary  fiber  0.5  mm  in 
diameter. 
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choice  of  magnification,  demagnification,  or  separation  to  discrete  chips.  In  addition,  a 
pre-patient  optic  could  be  employed  to  increase  the  available  intensity  in  a  fan  beam 
relative  to  simple  slot  collimation. 

Digital  X-ray  Detection 

In  almost  all  radiological  systems,  detection  of  x  rays  is  performed  by  the  use  of 
a  "phosphor"  screen  which  converts  absorbed  x  rays  into  visible  light  photons.  The 
visible  light  is  then  recorded  on  an  analog  medium,  such  as  film,  or  detected  by  a  digital 
detector.  The  phosphor  screen  is  used  because  most  detectors,  including  film,  are  not 
particularly  sensitive  to  x  rays.  The  use  of  a  screen-film  combination  reduces  radiation 
requirement  compared  to  direct  exposure  film  by  a  factor  of  50-1 00.7  This  increase  in 
sensitivity  occurs  at  the  expense  of  resolution.  The  modulation  transfer  function  of  a 
good  film-screen  combination  drops  to  less  than  0.1  at  a  frequency  of  15  Ip/mm.7  The 
choice  of  screen  thickness  is  a  trade-off  between  detected  quantum  efficiency,  which 
improves  with  increasing  thickness,  and  resolution,  which  is  degraded  by  light  blur  in  a 
thick  screen.8 

Efficient  direct  x-ray  detection  virtually  eliminates  the  tradeoff  between  spatial 
resolution  and  sensitivity  because  of  the  elimination  of  the  phosphor  screen.  Direct 
x-ray  detectors  can  provide  resolution  of  20  Ip/mm  or  better  with  nearly  100%  Detector 
Quantum  Efficiency  (DQE). 

An  especially  promising  direct  x-ray  detector  is  the  cadmium  zinc  telluride  (CZT) 
semiconductor  detector  recently  developed  by  Digirad,  formerly  Aurora  Technologies. 
The  DQE  of  a  .5  mm  thick  CZT  detector  is  essentially  unity  at  20  keV,  while  the 
resolution  is  to  first  order  independent  of  thickness.  CZT  detectors  are  similar  in  many 
respects  to  cadmium  telluride  (CdTe)  detectors  which  have  been  available 
commercially  for  many  years.  The  replacement  of  a  fraction  of  the  Cd  with  Zn  causes  a 
wider  bandgap  and  results  in  a  resistivity  increase  of  two  orders  of  magnitude.  This 
high  resistivity  is  an  important  factor  because  it  reduces  leakage  current,  a  significant 
source  of  performance  degrading  noise.  Leakage  current  noise  generally  limits  the  use 
of  CdTe  detectors  by  requiring  longer  integration  times  and  larger  pixel  sizes. 

Another  promising  direct  digital  detector  is  the  Charge  Injection  Device  (CID). 

Like  the  more  common  CCD  (Charge  couple  device)  CID  technology  allows  for  rapid 
array  imaging.  Unlike  the  CCD,  the  CID  is  very  radiation  resistant.  Also  unlike  CCDs, 
for  which  the  entire  array  must  be  read  out  and  simultaneously  cleared  as  a  unit,  CID 
pixels  can  be  read  individually,  repeatedly,  and  nondestructively,  leading  to  very  high 
dynamic  range.  CID  pixel  well  depths  are  also  deeper  than  conventional  CCDs,  leading 
to  a  substantial  reduction  in  blurring  from  high  energy  x-ray  photons. 
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Proposed  Statement  of  Work:  Status 

I.  Development  of  Capillary  Optic  Mammographic  System:  Months  1-36 

I.A.  Single  Capillary  Measurements:  Months  1-6:  Completed 

I.B.  Simulation  of  Design  Strategies:  Months  6-18:  Completed 

I.C.  Assembly  and  Testing  of  Prototype  Optic:  Months  18-24:  Completed 
I.D.  Design  of  Final  Mammographic  Optics  and  Detector  Unit:  Months  20-28: 
90%  Completed 

I. E.  Assembly  and  Testing  of  Final  Optic:  Months  28-36:  Delayed,  Ongoing 

II.  Development  of  Digital  Detector:  Months  1-24 

I  I.A.  Fabricate  CZT  Linear  Detector  Arrays:  Months  1-24:  Completed 

II.  B.  Develop  Interconnect  Methodology  (Wire  Bonding):  Months  2-8:  Completed 

II. C.  Assemble  and  Test  Proof  of  Principle  Detector/Multiplexer  Hybrid  :  Months 

2-12:  Completed 

II.  D.  Assemble  Additional  Hybrid  Arrays:  Months  12-24:  Changed 

III.  Mammographic  Measurements:  Months  18-48 

III.  A.  Measurement  of  Prototype  System:  Months  18-36:  Completed 

III.A.1  Design  and  Fabrication  of  Test  System:  Months  18-24:  Completed 

III.A.2.  Measure  Primary  Transmission  in  a  Mammographic  Geometry: 
Months  24-30:  Completed 

III.A.3  Measure  Transmission  of  Scattered  Radiation:  Months  24-30: 
Completed 

III.A.4.  Spectral  Measurements:  Months  24-30:  Completed 
III.A.5.  Investigate  Artifacts  due  to  Capillary  Structure:  Months  30-36: 
Ongoing 

III  .A.  6.  Measure  Contrast  Improvement  and  Resolution:  Months  30-36: 
Completed 

III.B.  Measurement  of  Final  Optics/Detector  System:  Months  28-48 

III.B.1.  Design  and  Fabrication  of  Measurement  System:  Months  28-36: 
Completed 

III.B.2.  Measure  Primary  Transmission  in  a  Mammographic  Geometry: 
Months  36-42:  Ongoing 

III.B.3.  Measure  Transmission  of  Scattered  Radiation:  Months  36-42: 
Ongoing 

III.B.4.  Spectral  Measurements:  Months  36-42:  Ongoing 
III.B. 5.  Investigate  Artifacts  due  to  Capillary  Structure:  Months  42-48: 
Ongoing 

III.B.6.  Measure  Contrast  Improvement  and  Resolution:  Months  42-48: 
Ongoing 

III.B.7.  Evaluate  Image  Quality  Using  RMI  Breast  Phantom:  Months  42- 
48:  Future 


t  t 
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Progress  to  Date 

Rather  than  follow  a  step  by  step  task  listing,  this  progress  report  has  been 
divided  into  subject  areas  to  provide  a  more  coherent  description  of  the  current  state  of 
the  art.  While  there  has  been  substantial  manufacturing  delays  in  the  production  of 
monolithic  optics,  the  available  optics  have  lent  themselves  to  a  thorough  systematic 
review  of  the  various  optics  geometries  proposed  for  this  project.  Representative 
measurements  have  been  made  on  a  large  number  of  optics  geometries,  rather  than 
relying  on  the  simple  simulations  initially  proposed  for  this  project.  Manufacturing 
issues  are  discussed  as  well.  A  decision  has  been  made  to  concentrate  on  two 
geometries,  a  collimating  optic  with  a  prepatient  monochromatizing  crystal,  and  a  “multi¬ 
taper”  post  patient  optic  assembled  of  three  smaller  tapers  to  produce  a  larger  area 
optic  without  the  technological  problems  of  a  single  large  area  glass  bundle. 
Manufacturing  delays  have  led  to  a  granted  no  cost  extension  of  this  project. 

The  proof-of-principle  linear  CZT  array  has  been  demonstrated,  but  an  alternate 
detector  path,  CID  technology  appears  equally  promising.  Studies  of  both  detectors  are 
discussed. 

Polycapillary  Optics 

Measurements  and  simulations  have  been  performed  on  several  geometries: 
single  fibers;  multifiber  collimating  optics;  multifiber  straight  antiscatter  optics;  and 
magnifying  and  demagniying  monolithic  optics. 

Single  Fibers 

Transmission 

An  important  step  in  optics  development  was  the  breakthrough  in  modeling 
capacity,  which  has  led  to  an  unprecedented  level  of  theoretical  understanding  of  the 
basic  properties  of  polycapillary  fibers,  especially  fiber  defects.  To  evaluate  the 
experimental  performance  of  polycapillary  fibers,  and  design  capillary  optics,  it  is 
necessary  to  be  able  to  predict  theoretical  behavior  for  complex  geometries.  The 
modeling  program  for  single  fibers  is  based  on  a  Monte  Carlo  simulation  of  simple 
geometrical  optics.  The  computational  speed  is  greatly  enhanced  by  a  reduction  to  two 
dimensions  by  projecting  the  trajectory  onto  the  local  fiber  cross-section.9  Reflectivities 
are  computed  from  standard  tables.10  Significant  recent  progress  has  been  made  in 
understanding  the  effect  of  capillary  profile  error,  waviness,  and  roughness  on  the 
transmission  spectra.11  This  is  extremely  important  in  providing  feedback  to  the 
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Type 

Description 

Outer  Diameter, 
mm 

Channel  Size, 
pm 

Open  Area 

Length, 

mm 

A 

Borosilicate 

0.5 

12 

65% 

105 

B 

Lead  glass 

0.5 

12 

52% 

95 

C 

Borosilicate 

0.75 

22 

50% 

136 

D 

Borosilicate 

4 

12 

55% 

130 

E 

Borosilicate 

0.3 

4-5 

55% 

105 

Table  1.  Description  of  polycapillary  fibers. 


■  exp,  Fiber  A  . sim,  Fiber  A 

A  exp,  Fiber  B 


manufacturing  process.  In  the  last 
year,  it  has  been  found  that  extremely 
good  fits  can  be  produced  with  only 
two  fitting  parameters  if  a  more 
physical  model  of  waviness  is 
employed.12 

The  experimental  arrangement 
for  single  capillary  measurements  is 
shown  in  Figure  3  and  Figure  4.  An 
optical  rail  affixed  to  an  optical  table 
carries  an  x-ray  source,  fiber  platform, 
and  x-ray  detector.  Each  can  be 
positioned  independently  in  three 

■  data 


Figure  5.  Measured  transmission  versus  energy  for 
polycapillary  fibers  listed  in  Table  L 

dimensions.  A  collimator  is  placed  before  the  fiber  and  any 
remaining  x-ray  leakage  around  the  fiber  is  eliminated  with 
metal  powder  or  filings.  The  measured  fibers  are  described 
in  Table  1. 

The  results  of  transmission  studies  as  a  function  of 
photon  energy  are  shown  in  Figure  5.13,14  All  of  the  fibers 
except  the  lead  glass  have  transmissions  at  20  keV  nearly 
equal  to  their  fractional  open  area  (the  fraction  of  the  cross 
section  of  the  capillary  which  is  open  space,  the  rest  being 
glass  walls).  This  transmission  corresponds  to  the  primary 
transmission  expected  for  a  linear  capillary  optic  employed 
as  an  antiscatter  grid. 


—  z=0 

—  z-Q.b  nm 


Figure  6.  Source  scans  at  10 
keV ! 
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■  data 
_z=0 
_ z=1  nm 


Source  position(mm) 


Figure  7.  Source  scans 
at  68  keV. 


The  transmission  falls  off  at  higher  energies.  Modeling  this 
transmission  spectra  has  lead  to  a  greater  understanding  of 
defects  in  polycapillary  optics. 

In  Figure  6,  simulations  with  or  without  roughness 
corrections  are  compared  with  the  experimental  data.  In  these 
measurements  the  source  is  scanned  transverse  to  the  fiber 
axis.  The  simulation  with  a  roughness  height  of  0.5  nm  fits  the 
experimental  data  quite  well.  It  is  definitely  over-corrected  when 
the  roughness  is  1.0  nm.  The  same  simulations,  shown  in 
Figure  7,  are  also  carried  out  at  68  keV,  where  the  width  of  the 
curve  is  narrower  than  that  at  1 0  keV  because  of  the  smaller 
critical  angle.  As  a  result,  photons  also  experience  fewer 
bounces  on  average.  Simulations  with  roughnesses  as  large  as 
1 .0  and  2.0  nm  still  could  not  fit  the  data.  Knowing  that  1 .0  or 
2.0  nm  roughness  is  definitely  too  large  at  10  keV,  we  can 
determine  that  the  roughness  correction  by  itself  is  not  sufficient 
at  high  energy  to  reproduce  source  scan  measurements.  Other 
effects  need  to  be  considered.  These  are  bending  and 
waviness. 


A  slight  bending  of  the  capillary  can  dramatically  reduce  the  transmission  of  high 
energy  photons  because  of  the  small  critical  angle.  A  comparison  between 
experimental  data  and  simulations  with  different  bending  is  shown  in  Figure  8.  The 
figure  shows  that  the  simulations  with  bending  alone  do  not  fit  the  experimental  data 
well,  which  indicates  that  bending  is  not  the  only  factor  which  causes  the  high  energy 


0s* 


Figure  8.  Transmission  spectra  of  Fiber  3  simulated  waviness  only  compared  with  the  experimental  data, 
with  different  bending  curvature  alone  and  compared 
with  experimental  data 
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Source  position(mm)  Source  position(mm)  Source  positjon(rrm)  Source  position(mm) 


Figure  10.  Simulated  source  scan  curves  compared  with  experimental  data  at  four  different  photon 
energies.  Parameters  are:  R  =  125  m,  AOmax=  0.35  mrad,  roughness  height-0.5  nm. 


transmission  to  drop.  However,  from  Figure  8,  we  can  see  that  the  range  of  the  bending 
radius  must  be  larger  than  100m  to  give  the  observed  transmission  at  the  highest 
energy(80KeV). 

Capillary  surface 


oscillations  with  wavelengths 
shorter  than  the  capillary 
length  and  longer  than  the 
wavelength  of  the  roughness 
are  called  waviness.  The 
detailed  shape  of  waviness  is 
unknown.  Its  average  effect 
can  be  considered  as  a 
random  tilt  of  the  glass  wall,  so 
that  the  grazing  angle  of  the 
photon  is  changed  by  a 
random  amount,  SO,  after 
every  bounce.  SO  is  a  random 
number  between  -A6may  and 
^9  max  if  9  >  AO  max.  The 


maximum  random  tilt  angle 
A0max  is  an  adjustable 


Photon  energy(Ke\/) 


parameter  which  depends  on  Figure  11.  Simulated  transmission  spectra  with  different 


the  waviness  of  the 


bending  and  waviness  compared  with  the  experimental  data  in 


polycapillary  fiber  To  keep  O'  search  for  the  best  fitting  of  Fiber  C. 
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Figure  12.  Scheme  of  an  x-ray  beam  interacts  with  a 
random  tilt  surface. 


experimental  data. 


positive,  SO  is  taken  to  be  a 
random  number  between  -0 
and  A0max  when  0  <  A6max . 

Since  a  photon  with  a  incident 
angle  smaller  than  A0max  has  a 
larger  chance  to  experience  an 
angle  increase  than  an  angle 
decrease,  this  is  physically 
reasonable.  In  Figure  9, 
simulations  with  waviness 
corrections  with  A0max  set  at  1 
mrad  and  2  mrad,  which  is 
comparable  to  the  critical  angle, 
are  compared  with  the 

This  figure  shows  that  simulations  with  waviness  alone  do  not  fit  the 

experimental  data.  This  is  because  the 
waviness  correction  changes  the 
reflected  angle,  not  the  profile.  In  fact 
the  capillary  is  still  considered  to  be 
straight,  so  those  photons  which  have 
few  reflections  will  not  be  significantly 
effected  by  waviness. 

Finally  the  waviness  and 
bending  are  combined  by  increasing 
the  bending  radius  R,  roughly 
determined  in  Figure  8,  and  adding  a 
waviness  parameter,  AGmax.  Several 
trials  are  shown  in  Figure  1 1 .  Sim2  has 
too  much  bending  and  not  enough 
waviness;  sim3  has  too  much  waviness 
and  not  enough  bending;  siml  is  the 
best  fit.  Roughness  is  also  included  in 
those  simulations.  The  source  scan 

simulation  with  the  three  fixed  parameters  are  plotted  along  with  the  experimental  data 
in  Figure  10  for  four  more  photon  energies.  They  all  fit  quite  well.  This  model,  hereafter 
labeled  Ml,  assumed  a  uniform  distribution  of  tilt  angles,  which  is  not  quite  physical.  In 
more  recent  work,  discussed  as  Model  M2,  it  is  assumed  that  these  tilt  angles,  <j>,  are 
normally  distributed  in  the  range  (-tc/2,  n/2)  with  the  mean  value  equal  to  zero.  For  high 
quality  optics,  the  standard  deviation  of  this  normal  distribution,  a,  is  much  smaller  than 
the  critical  angle,  Gc.  The  probability  distribution  of  tilt  angles,  4>,  is 


Figure  13.  Three  surfaces,  OA1,  OA2  and  OA3,  with 
different  tilt  angles,  <j)l,  </>2  and  <j>3,  respectively,  from  the 
nominal  surface  OO  ’. 


G(<t>)  = 


1 


W 27 


r 

la2 


(1) 
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In  this  work,  Model  M2,  consideration  was  taken  of  the  fact  that  the  surface  tilt  angle  will 
affect  the  probability  of  x-ray  impact  on  that  surface.  Taken  to  extremes,  a  surface 
region  perpendicular  to  the  beam  is  much  more  likely  to  intercept  the  beam  than  a 
surface  region  parallel  to  the  beam.  Figure  13  displays  three  surfaces,  OA^  OA2  and 
OA3,  with  different  tilt  angles,  <J>1t  <|)2  and  <J>3,  respectively,  from  the  nominal  surface  00’. 
The  projections  onto  the  nominal  surface  for  the  three  surfaces  are  equal,  OAj *cos<t>.,  = 
OA2*cos<j>2  =  OA3*cos<t>3=  00’.  The  probability  of  incidence  of  a  parallel  x-ray  beam  with 
incident  angle  0  (with  respect  to  the  nominal  surface  00’)  hitting  these  tilted  surfaces  is 
given  by  their  corresponding  perpendicular  length  AjBj  (for  j  =  1,  2  or  3) 

P,  ccA.B,  =OA  ■  sin(<9  +  <f>j ) ■  sin(6l  +  4 . ),  ^ 

COS(f>j 

We  call  this  the  tilt-corrected  probability  distribution.  The  complete  description  of  this 
distribution  is 


F  sin(  0  + , 
cos( 


where  0  is  the  incident  angle,  <|>  is  the  tilt  angle  and  F  is  a  normalization  constant. 

Combining  the  normal  distribution  G(<t>)  with  the  tilt  correction  H(<|>)  gives  <|>  for  a 
certain  incident  angle  0  as  follows: 

kJ$'2&+£>'  _e<ip<£ 

J{(j))  =  \  COS0  2, 

0, 

where  0  is  the  incident  angle,  <|>  is  the  tilt  angle  and  K  is  a  normalization  constant. 
Noting  that 

5^^+corf-ta*  <5) 

COS^ 

I 

and  that  both  the  incident  angle,  0,  and  the  tilt  angle,  <f>,  are  very  much  less  than  40 
mrad,  we  use  the  approximation 

f  dtL 

Ce  *  {6  +  (j>), 

P{*)  =  2  , 

0, 

where  the  normalization  constant ,  C,  is 


yJlnOcr  +  2cr 2  ’ 


(7) 
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Figure  14.  P(</>)  and  J(</>)  versus  tilt 
angle  <j>  when  incident  angle  0  = 
0.009rad,  standard  deviation  o- 
0.0004rad.  <f>  is  in  rad. 


Figure  15.  Simulations  of 
transmission  spectra  for  fiber  C  with 
only  waviness  compared  with  the 
experimental  data.  This  figure  shows 
the  effects  of  waviness.  It  also  shows 
simulations  using  waviness  alone  do 
not  fit  the  data.  The  simulations  do 
not  include  roughness  or  bending. 


to  speed  the  calculation.  P(<f>)  is  called  the  tilt-corrected 
normal  distribution.  In  model  M2,  a  is  the  variable  that 
describes  the  amount  of  waviness.  An  example  of  J(<j>)  and 
P( ([>)  is  shown  in  Figure  14.  The  result  of  the  waviness 
calculation  is  shown  in 

Figure  15.  The  model  M2  shows  a  much  larger 
drop  at  high  photon  energies  for  lower  values  of 
the  waviness  than  does  model  Ml.  A 
comparison  of  the  simulation  with  bending  and 
waviness  to  the  experimental  data  is  shown  in 
Figure  16. 

Using  a  similar  technique,  the  simulation 
was  compared  to  the  fibers  listed  in  Table  1 , 
using  the  parameters  listed  in  Table  2.  The 
results  are  plotted  in  Figure  17  and  Figure  18.  In 
Figure  1 7,  the  transmission  for  fiber  A  shows  a 
rapid  drop  for  energies  above  30  keV.  Although 
fiber  D  has  lower  fractional  open  area  than  fiber 
A,  its  transmission  exceeds  that  of  fiber  A  at 
energies  above  30  keV.  This  is  because  fiber  A 


Figure  16.  Simulated  transmission  spectra  using 
model  M2  with  fixed  waviness  ((7=0.2  mrad)  and 
different  bending  (bisection  method),  compared 
with  the  experimental  data.  The  simulation  curves 
shown  are  just  a  few  representative  cases  in  the 
bisection  process.  For  M2,  the  best-fit  parameters 
are  a  =  0.2  mrad  and  R  -  225  m;  for  Ml,  they  are  z 
=  0.5  nm,  s  =  6  pm,  AQmax  =  0.35  mrad  and  R  = 
125  m. 
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parameters  compared  with  the  experimental 
data.  All  parameters  related  to  this  figure  are 
listed  in  Table  1. 

transmission  about  the  same  as  for  fiber  C. 


is  thin  (0.5  mm  in  outer  diameter)  and 
flexible,  therefore  difficult  to  keep 
straight  in  the  measurement  apparatus. 
Model  M2,  which  can  vary  only 
bending  and  waviness,  requires  a 
much  sharper  bend  for  fiber  A  than  for 
fiber  D.  This  bending  is  more 
significant  at  high  energies,  where  the 
critical  angles  are  smaller. 

Transmission  curves  for  fiber  C  and  D 
are  similar  in  shape  and  are  nearly  flat 
up  to  60  keV,  as  shown  in  Figure  18. 
However,  the  smaller  channel  size  for 
fiber  D  results  in  the  simulation  being 
less  sensitive  to  bending.  Therefore  a 
larger  bending  curvature  (smaller 
radius)  is  required  for  fiber  D.  The 
waviness  correction  for  fiber  D  must 
therefore  be  smaller  to  keep  the 


Figure  17.  Simulations  of  transmission  spectra 
of  fiber  A  and  fiber  D  with  their  best-fit 


Fiber  E  is  the  thinnest  fiber.  Table  2  shows  that  it  has  the  largest  bending 
correction  because  of  its  flexibility.  Its  transmission  curve  is  flat  up  to  40  keV.  This  is 
due  to  its  very  small  channel  size,  as  discussed  in  the  comparison  of  fiber  C  and  fiber 

D.  However,  if  the  channel  size  is 
too  small,  as  in  fiber  E,  it  also 


Fiber  C:  X  data,-  -  Ml :  best  fit*  A  M2xf=0.2,R=225 

70 1  Fiber  D:  X  data, -  >-  •  Ml:  best  M2:o=0.125,  R=139 

‘  Fiber  E:  -+ . data,— ▼—  Ml:  Best  frtrO—  M2:cr=0.09,  R=3 1 

60  . 


Photon  energy  (KeV) 


results  in  more  reflections  being 
needed  for  a  photon  to  traverse 
the  fiber  and  may  have  introduced 
other  defects  such  as  blocked 
channels.  This  is  why  the 
transmission  is  only  40%  for  the 
energies  below  40  keV  although 
the  open  area  is  around  55%. 

The  new  model  yields  as  high  a 
quality  fit  with  two  fewer  fitting 
parameters  than  used  in  model 
Ml  in  Figure  5.  The  high 
transmission  and  the  simulation 


Figure  1 8.  Simulations  of  transmission  spectra  of  results  show  that  the  quality  of  the 
fiber  C,  fiber  D  and  fiber  E  with  their  best-fit  capillary  fibers  is  quite  good.  The 

parameters  compared  with  the  experimental  data.  bending  radius  is  above  130 

All  parameters  related  to  this  figure  are  listed  in  meters  for  type  3  and  type  4 
Table  2.  capillaries.  It  is  hopeful  that  we 
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Fiber  Description 

Model  Ml 

Model  M2 

Fiber 

# 

Type 

Outer 

Diameter 

(mm) 

Channel 

Size 

(um) 

Area 

Length 

(mm) 

Z 

nm 

s 

pm 

R 

m 

E  (0 

<x>  i= 

<1  E 

R 

m 

CT 

Mrad 

A 

1 

0.5 

12 

65% 

105 

0.7 

6 

105 

0.4 

60 

0.225 

C 

3 

0.75 

22 

50% 

136 

0.5 

6 

125 

0.35 

225 

0.2 

D 

4 

4 

12 

55% 

130 

0.8 

6 

110 

0.285 

139 

0.125 

E 

5 

0.3 

MEM 

55% 

105 

0.7 

6 

28 

0.2 

31 

0.09 

F 

4 

4 

12 

55% 

130 

0.8 

6 

90 

0.45 

90 

0.18 

Table  2  Parameters  for  best-fit  simulations.  R  is  the  bending  radius.  For  Ml,  A0maxis  the 
amount  of  waviness,  z  is  the  roughness  height  and  s  is  the  roughness  correlation  length. 
For  M2,  cris  the  standard  deviation  of  the  waviness. 


can  further  improve  the  high-energy  transmission  performance  of  polycapillary  fibers  by 
decreasing  to  an  optimum  channel  size,  and  making  them  more  rigid. 


Fiber  B  in  Figure  5  has  poor  transmission;  it  is  a  lead  glass  fiber,  with  high 
surface  roughness  and  waviness.  More  recently,  higher  quality  lead  glass  fibers  have 
been  produced.  The  transmission  as  a  function  of  photon  energy  for  these  fibers  is 
shown  in  Figure  19.  These  leaded  fibers  have  higher  absorption  of  high  angle  photons 
than  the  lower  density  borosilicate 
glass  fibers,  and  therefore  would 

provide  even  better  scatter  rejection.  x  F 


Transmission  was  also 
measured  of  deliberately  curved  fibers, 
which  can  be  used  to  deflect  x-ray 
beams  over  larger  angles.  The  results 
are  shown  in  Figure  20.  Simulated 
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Figure  20.  Transmission  efficiency  of 350  pm  OD 
polycapillary  fiber  with  6  pm  channel  diameter .  The 
fiber  length  is  25  cm.  Simulated  values  are  shown  for 
20 ,  30,  40,  50  keVxrays.  The  parameters  are 
curvature,  R;  waviness ,  w;  and  roughness,  z. 


Figure  19.  Recently  manufactured  lead  glass 
fibers. 
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Figure  21.  Simulated  transmission  as  a  function  of 
bending  radius  at  50  keV,  with  no  waviness  or  roughness. 


transmission  for  a  single  fiber  with 
various  channel  sizes  is  shown  in 
Figure  21. 

Absorption  Measurement 


To  achieve  total  scatter 
rejection,  the  glass  wall  of  the  optics 
must  be  thick  enough  to  absorb  the 
scattered  x  rays  from  various 
directions  while  allowing  the  majority 
of  the  primary  x  rays  to  pass.  The 
absorption  of  the  polycapillary  fiber 
bundle  was  measured  by  tilting  the 
optic  to  an  angle  which  is  larger  than 


the  critical  angle. 

The  setup  includes  a  source,  a  pinhole  with  size  of  200  pm,  an  aluminum  plate 
which  is  used  to  hold  the  capillary  bundle,  and  a  Germanium  detector.  The  source  has 
a  maximum  operate  voltage  of  100  kV.  The  source  and  aluminum  plate  are  attached  on 
linear  stages  which  are  controlled  by  linear  actuators.  Both  the  source  and  the  optics 
fiber  can  move  freely  in  two  perpendicular  directions  transverse  to  the  x-ray  beam.  An 
additional  rotary  stage  is  used  to  rotate  the  aluminum  plate,  so  that  the  optic  can  also 
be  rotated  in  horizontal  direction. 

The  experiment  starts  with  the  alignment  between  the  source,  pinhole,  and 
detector,  by  adjusting  the  source  and  fiber  positions.  After  alignment,  the  optic  was 
tilted  to  an  angle  larger  than  the  largest  critical  angle  for  the  measured  energy  range, 
then  translated  out  of  the  x-ray  beam.  While  the  capillary  optics  was  moved  into  the 
x-ray  beam  step  by  step,  an  x-ray  spectrum  was  taken  by  the  germanium  detector  after 

the  optics  at  each  step.  The  acquire  time 
for  each  step  was  30  seconds.  The 
transmission  at  every  step  is  the  ratio 


o 
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Figure  22.  Absorption  measurement 


Figure  23.  High  angle  transmission  through 
capillary  walls  as  a  function  of photon 
energy. 
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between  the  counts  without  optics  and  the  counts  with  optics  at  that  step.  Transmission 
at  different  energies  as  a  function  of  optic  position  is  plotted  in  Figure  22  along  with  the 
theoretically  calculated  values. 

Large  angle  transmission  can  be  calculated  as 

J  _  e~(l (2) 
s 


where  p  is  the  mass  attenuation  of  the  material  which  can  be  calculated  from  tabulated 
values21  according  to  the  mass  composition,  p  is  the  density  of  the  glass,  f  is  the  open 
area  of  the  capillary  bundle,  and  x  is  the  distance  that  the  x-ray  beam  travels  across  the 
fiber,  x  can  be  calculated  from  capillary  movement  d  and  capillary  tilting  angle, 
x  =  d/ a ,  /?) 


Lead  glass 


Figure  25.  Calculated  scatter  transmission  for  lead 
glass  optics  of  various  lengths. 


where  a  is  the  tilting  angle  in  rad. 
Theoretical  calculated  numbers  from 
equation  (2-3)  are  compared  with  the 
experimental  result  in  Figure  22.  The 
theoretical  values  fit  well  with  the 
experimental  data.  Measured  and 
theoretical  transmission  of  a 
polycapillary  fiber  as  a  function  of 
photon  energy  is  shown  in  Figure  23. 
This  "high  angle  transmission" 
corresponds  to  scatter  transmission  for 
the  capillary  optics  used  as  an 
antiscatter  grid,  and  is  less  than  1%  at 
20  keV.  The  solid  lines  are  theoretical 
calculations,  which  agree  quite  well 
with  the  data.  Lead  glass  optics  would 
have  even  better  scatter  rejection,  as 
shown  in  Figure  25. 15 


Collimating  Optics 

8  keV  Prototype 

Transmission  and  Uniformity 

Because  of  the  technological  difficulties  in 
producing  large  area  monolithic  optics,  the 
fastest  route  to  a  large  scale  antiscatter 
system  for  contrast  measurements  is  a 
straight  optic. 


Figure  24.  Multifiber  collimating  lens. 
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The  low  angular  acceptance  of  the  capillary  channels  requires  that  such  an  optic 
be  used  with  a  collimated  beam,  as  shown  in  Figure  26.  A  collimating  lens,  similar  to 
that  pictured  in  Figure  24,  and  a  straight  fiber  bundle  of  3  cm  diameter  have  been 
obtained.  The  collimating  lens  has  an  output  diameter  of  3  cm,  a  focal  distance  of  150 
mm,  and  collects  from  an  angle  of  8°.  It  was  designed  for  8  keV,  where  the 
transmission  is  30%,  but  has  a  transmission  of  14%  at  20  keV.  The  drop  off  is  due  to 
the  reduction  in  critical  angle  for  reflection  with  x-ray  energy,  as  stated  in  equation  (1). 
Because  this  optic  would  be  used  prepatient,  low  transmission  would  not  increase 
dose,  only  increase  tube  loading.  More  importantly,  as  will  be  demonstrated  in  the  next 
section,  a  lens  actually  designed  for  20  keV  would  have  a  higher  transmission  at  that 
energy.  The  primary  purpose  of  studying  this  lens  is  to  develop  experience  with  large 
scale  polycapillary 
optics  and  refine  the 
modeling  to  increase 
the  confidence  in  the 
simulations.  The 
simulated  and 
measured  transmission 
of  the  lens  as  a  function 
of  source  to  lens 
distance  is  shown  in 
Figure  27.  The 
simulation,  which 
models  the  complex 
geometries  of  complete 
lenses,  was  developed 
by  X-ray  Optical 
Systems,  Inc.  Some 
delay  in  the 
commencement  of  the 
simulation  task  was 
created  by  the  necessity 
of  assuring  the 
confidentiality  of  the 
proprietary  code  while 
providing  access  to 
student  researchers. 

Access  has  been 
provided,  and 
comparisons  of 
measurements  and 
simulations  of  a 

Collimating  lens  have  figure  27.  Variation  in  measured  (box)  and  simulated  (x)  20  keV 
given  goo  resu  S.  e  transmission  as  the  source  is  moved  towards  the  lens  along  the  fiber 

axis. 


Figure  26.  Set-up  for  Anti-scatter  measurements. 
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measured  and 
simulated 
transmission  as 
a  function  of 
photon  energy 
for  the  lens 
(with  the 
central  part 
blocked,  as 
Figure  28.  X-ray  photograph  of  shown  in 
output  of  collimating  optic.  Figure  29),  is 

shown  in 


Figure  29.  X-ray  photograph  of 
output  of8keV collimator.  A 
rectangular  lead  sheet  blocks  part  of 
the  output. 


Figure  31. 

More  serious  for  imaging 
experiments  than  the  low 
transmission  of  a  lens  designed 
for  lower  energies,  is  the  non¬ 
uniformity  of  the  field.  The  outer 
fibers  are  bent  too  much  to  permit 
good  transmission  at  these 
energies.  The  fall  off  in 
transmission  can  be  seen  in  the 
photograph  shown  in  Figure  28, 

and  in  the  simulation  of  Figure  30.  Figure  30.  Simulated  output,  at  20  keV,  of  collimator 
Contrast  enhancement  designed  far  8  keV. 


measurements,  which 
are  ongoing,  must  take 
these  nonuniformities 
into  account.  Figure  29 
shows  the  output  of  the 
lens  after  a  center 
section  was  blocked. 

This  section  contains 
fibers  which  were 
removed  for  study  of 
radiation  exposure 
effects  after  extensive 
exposure  to  intense 
synchrotron  radiation.16,17 
The  radiation  effects 
were  found  to  be 
manageable  and  are  not 
expected  to  be 
significant  for 


I . T . I . I . . . I . T . 1 . 1 

14  15  16  17  18  19  20  21 

Energy  (KeV) 


Figure  31.  Transmission  of  8  keV  collimator  as  a  function  of photon 
energy.  The  simulation  assumes  a  fiber  open  area  of  65%,  the  fibers 
packed  to  fill  50%  of  the  lens  input  area,  and  ‘A  of  the  lens  blocked  by 
the  lead  shield.. 
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Xads,  rrm  source  scan,  as  a  function  of  source  motion 

along  the  lens  axis,  for  collimating  lens, 


Figure  33.  Transmission  of  a  multifiber  collimating  measured  at  20  keV. 
lens  as  the  source  is  moved  transverse  to  the  lens  axis. 

Circles  are  measured  data,  and  boxes  simulation 
results. 


mammography. 

The  aligned  transmission  measurements  of  Figure  27  and  Figure  31  are 
important  to  determine  the  collection  efficiency  of  the  lens.  Variation  in  transmission  as 
a  function  of  entrance  angle  is  important  for  estimating  the  scatter  rejection  properties 
of  other  optics  which  might  be  used  post  patient.  This  variation  was  measured  by  first 
aligning  the  source  with  the  lens  and  then  moving  the  source  in  the  plane  transverse  to 
the  capillary  axis.  The  transmission  efficiency  decreases  when  the  x-ray  source  is 
moved  off-axis  because  the  angle  between  the  x  ray  and  the  reflection  surface  of  the 
capillary  increases.  Because  the  reflectivity  decreases  sharply  with  angle  this  produces 
greater  loss  at  each  reflection.  In  addition,  x  rays  traveling  through  the  capillary  at  a 
larger  angles  undergo  more  reflections.  A  comparison  of  measured  and  simulated 
data,  with  good  agreement,  is  shown  in  Figure  33.  Variation  in  the  width  of  the  scan  as 
a  function  of  the  lens  to  source  distance  is  also  used 


as  the  primary  method  of  determining  both  the  focal 
distance  of  the  lens,  and  the  depth  of  field.  The 
measured  variation  is  plotted  in  Figure  32.  The  depth 
of  field  can  also  be  seen  in  Figure  27,  which  shows  the 
measured  and  simulated  transmission  as  a  function  of 
source  to  lens  distance.  The  depth  of  field  is  about  1 1 
mm. 

Divergence  and  Diffraction  Gain 

Measurements  of  the  exit  divergence  of  the 
collimating  optic  have  been  performed  at  8  keV  and  at 
20  keV.  Exit  divergence  is  important  if  a  capillary  optic 
is  to  be  used  as  a  collector  fore  slit  before  the  patient. 


Figure  34.  Image  of  lens 
shown  in  figure  30,  taken 
30  cm  after  the  lens. 
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In  this  case,  the  spatial 
resolution  will  depend  on  the 
angular  spread  of  x-rays  at  the 
exit  end  of  the  polycapillary 
fibers.  Large  angular  divergence 
is  equivalent  to  a  large  focal  spot 
size  in  a  conventional  system. 
The  divergence,  measured  by 
rotating  a  silicon  wafer  in  the 
beam  and  measuring  the  width 
of  the  diffraction  curve  is  shown 
in  Figure  35.  The  full  width  of 
the  divergence  measured  at  20 
keV  is  2.56  mrad,  slightly  greater 
than  the  critical  angle  for  total 
external  reflection.  The  effect  of 
the  divergence  is  shown  in  figure 
34,  in  which  the  image  of  the  fiber  bundle  walls  have  been  blurred  by  placing  the  film  30 
cm  behind  the  lens. 

The  ability  of  the  optic  to  redirect  the  collected  x  rays  into  a  beam  of  relatively 
low  divergence  means  that  the  reflection  efficiency  off  a  monochromatic  crystal  is 
greatly  enhanced.  While  the  low  power  test  source  was  not  intense  enough  to  permit 
imaging  after  monochromatizing,  measurements  were  made  of  crystal  reflectivity.  A 
silicon  (100)  crystal  was  found  to  reflect  approximately  .0047%  of  the  beam  from  the 
molybdenum  source  after  collimation.  This  is  in  reasonable  agreement  with  a 
calculation18  of  expected  reflectivity  for  a  crystal  with  0.14  mrad  mosaicity,  which  is 
reasonable  for  the  crystal  used.  The  gain,  relative  to  no  collimation,  is  on  the  order  of 
lOOx. 


20  KEV 
PROTOTYPE 

A  second 
collimating  lens 
designed  for  higher 
energies  was  also 
tested.  This  lens  has 
a  focal  length,  1m, 
which  was  longer  than 
ideal  for 

mammography,  but  a 
much  larger 
transmission,  in 
excess  of  30%  at  20 
keV,  as  shown  in 
Figure  36. 


Figure  35.  Output  divergence  of  the  8  keV  collimating  lens, 
measured  at  20  keV. 


Figure  36.  Simulated  and  measured  transmission  of  the  20  keV  prototype. 


Direct  Digital  Mammography  Using  Capillary  Optics 


24 


■  Experiment 

•V1  •  Double  Bending  Simulation 
O  Single  Bending  Sim  with  same  parameters 
X  Single  Bending  Best  Fit 


Figure  37.  X-ray  micrograph  of  20 
keV  collimator,  exposed  5  min  at  60 
keV  and  0.4  mA  with  a  W  source  and 
Polaroid film. 


Figure  40.  Transmission  efficiency  of  a  central  (nominally 

straight)  polycapillary  fiber  in  the  20  keV prototype  *  he  drop-off  at  higher  photon 

collimator.  The  fiber  OD  is  350  pm  and  the  channel  energies  is  due  to  some  curvature  of 

diameters  are  10  pm.  The  fiber  length  is  14  cm.  Simulated  eVen  the  Central  fiberS-  35  eXPlored  in  * 

values  are  shown  for  a  single  bend  and for  a  compound  be  Simulation  of  Figure  40.  The 

(double)  bend.  All  of  the  simulation  results  assume  a  output  field  W3S  also  much  more 

roughness  z  =  0.3  nm  and  waviness,  w  =  0.2  mrad.  uniform  than  for  the  previous 

collimator,  as  shown  by  the  picture  in 

Figure  37  and  the  measurement  of  Figure  39.  At  lower  exposure,  in  Figure  38,  it  can  be 
seen  that  there  were  some  defects  in  the  right  half  of  the  lens.  This  was  due  to  some 
fiber  misalignment,  which  was  reflected  in  the  divergence  measurements.  This  lens 
was  a  first  test  of  anew  screen  technology,  which  has  since  been  improved. 


0.4500  -  0.5000 
0.4000  -  0.4500 
0.3500  -  0.4000 
0.3000  -  0.3500 
0.2500  -  0.3000 
0.2000  -  0.2500 
0.1500  -  0.2000 
0.1000  ~  0.1500 
0.0500  -  0.1000 
0  -  0.0500 


Figure  38.  26—  ray  photograph  of  20 
keV  collimator,  exposed  5  min  at 
35  keV  and  0.25  mA.. 


Figure  39.  Output  intensity  distributions  for  the  prototype  TYPE  I 
collector  for  different  energies. 


Direct  Digital  Mammography  Using  Capillary  Optics 


25 


Preliminary  imaging 
measurements  have  been 
performed  with  the  optic.  Figure  42 
and  Figure  41  are  radiographs  of  a 
5  cm  thick  Lucite  phantom.  Figure 
42  was  taken  with  the  film  5  cm 
from  the  phantom.  Figure  41  was 
taken  with  the  film  95  cm  from  the 
phantom.  There  is  no 
magnification  for  this  parallel  beam 


Figure  42.  Polaroid film 
radiograph  of  Lucite  phantom 
placed  5  cm  from  the  phantom. 
Upper  holes  are  6  mm  in 
diameter,  second  row  is  4  mm, 
third,  2  mm,  and fourth,  1  mm. 
Holes  on  the  right  are  3.5  cm 
deep,  middle  holes  are  3  cm 
deep  and  holes  on  left  are  2 
cm  deep. 


Figure  41.  Polaroid  film  geometry.  However,  some  blurring 

radiograph  95  cm  from  the  due  to  fiber  misalignment  is 

phantom.  observable.  Given  the  long 

distance,  known  fiber  misalignment 
and  very  low  exposure  the  image  is  reasonably  clear.  The 
images  were  taken  with  normal,  x-ray  insensitive  Polaroid 
film  and  a  very  low  power,  0.20  mA  test  source.  In  both 
cases  the  phantom  was  placed  close  to  the  output  of  the 
lens. 


2 

t- 


Anti-Scatter  Optics 
Straight 
MultiFiber 
Optic 

Rather  than 


simply  using  long  air 
gaps  to  eliminate 
scatter,  as  in  Figure  41, 
a  straight  polycapillary 
antiscatter  optic  could 
be  employed  ,  as 
shown  in  Figure  26. 
Such  an  anti-scatter 


optic,  which  was 
essentially  a  3  cm 


10  15  20  25  diameter  bundle  of 


Energy  (KeV)  straight  polycapillary 

fibers,  was  aligned  to 
the  output  of  the  8  keV 

Figure  43.  Simulation  of  straight  antiscatter  lens,  compared  to  measured  collimating  lens.  The 
data.  Simulation  uses  1  nm  roughness,  0.45  mrad  waviness,  and  26%  anti  -  scatter  optic  was 

overall  open  area.  1  50  mm  jn  |ength  wjth 

55|im  channels.  The 
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Figure  45.  2-D  source  scan  of  early  prototype 


optic  (input  diameter  =  2  mm,  output  diameter  =  3 
mm,  transmission  =  5  %  at  20  keV,  focal  distance 
22  cm)  This  lens  transmits  poorly  because  only  the 
central  capillary  channels  transmit  while  the  outer 
capillaries  do  not  transmit  at  all.  Only  about  33 
%  of  the  optic  cross  section  transmits  at  5%  or 
greater. 


transmission  spectrum  of  this  optic  is 


source  lens  distance(mm) 

Figure  44.  Source  scan  of  tapered  lens. 


shown  in  Figure  43.  Simulations  show  that  the  apparently  straight  fibers  have 
significant  curvature,  with  radius  of  curvature  of  about  25  m. 


Magnifying  Linear  Monolithic  Optics 

Linear  monolithic  optics  taper  as  a  single  glass  unit  from  input  to  output  as 
shown  in  Figure  2  and  Figure  46.  These  optics  are  relatively  difficult  to  manufacture. 
Many  early  prototypes  have  poor  transmission,  as  shown  in  Figure  45.  To  help 

understand  the 


Scanning  Geometry 


detector 


Figure  46.  Small  optic  scanning  geometry  for  contrast  and  imaging  measurements. 


defect  problem 
an  extensive 
simulation  and 
measurement 
study  was 
performed  on 
one  of  the  better 
prototypes. 

The  focal 
distance  of  a 
polycapillary 
optic  can  be 
determined  by 
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normalized 
transmission, 
exp  data 


Source  lens  distance  (mm) 


Figure  47 '  Transmission  of  tapered  lens. 


Transmission  spectrum  for  the  whole  lens 

)  ♦  Source  lens 

.  T  *  • — distance:  455  mm. 

3  t  4  \  - Simulated 

3  +  x\  transmission. 


finding  the  smallest  scan 

121  I  angle  or  the  maximum 

i  -  # — v  .  .  transmission.  Source  scan 

0.8-  g  ^X.  simulated,0"  curves  in  both  x  (horizontal) 

06  fj  \  normalized  and  y  (vertical)  directions 

of  ^  were  recorded  after  the 

0  4  source  and  the  lens  were  in 

0.2-/  ,.  .  alignment.  FWHM  of  the  scan 

o  J _ I _ I _ I  transmission,  curves  in  both  directions  are 

250  450  650  850  exP  data  calculated  with  Gaussian 

Source  lens  distance  (mm)  fitting,  and  the  average 

FWHM  is  taken  as  the  scan 
width  at  that  position.  Scan 

Figure  47.  Transmission  of  tapered  lens.  ang|e  js  defined  as  the  result 

of  the  average  scan  width 

Transmission  spectrum  for  the  whole  lens  divided  by  the  SOurce-lens  distance.  In 
40  T  ♦  source  lens  Figure  44,  scan  angle  is  plotted  as  a 

35  --  *  *XX. simulated455  mm'  function  of  source  lens  distance.  Since 

g  30  -  X.  transmission.  the  theoretical  curve  does  not  have  as 

§  25  -  ♦  \  many  data  points  as  the  experimental 

|  20  •  #  curve,  it  is  not  as  smooth  as  the 

|  15  "  ♦  X.  experimental  curve.  This  method  is  not 

2  10  "  ♦  ♦  ♦  *X^^  effected  by  the  source  instability,  and 

H  ®  ^  hence  is  very  useful  when  the 

0  0  20  40  60  80  experiment  needs  a  long  time,  and  the 

Photon  energy  (keV)  source  is  not  stable.  However  the  curve 

c,.  AO  „  , ,  is  quite  flat  near  the  focal  distance,  as 

expected  by  simulation.  The  scan 

widths  have  to  be  measured  very  carefully  to  determine  the  focal  distance.  The  focal 
distance  determined  by  minimizing  the  scan  angle  is  455  mm  for  this  lens. 

Transmission  and  photon  counts  are  also  plotted  versus  source  lens  distance  in  Figure 
47  along  with  the  simulation  values.  The  maximum  point  for  transmission  is  around  450 
mm.  The  second  method  can  be  affected  by  the  source  instability. 

Transmission  was  measured  in  a  wide  energy  range  for  the  whole  lens  as  well 
as  the  central  part  of  the  lens.  The 
transmission  spectrum  for  the  whole  _  70T  □ 

lens  is  plotted  in  Figure  48.  The  S  60”  *  pinhole  size:  i 

figure  shows  that  the  transmission  |  4*dD^X  mm 

drops  quickly  for  photons  with  energy  |  30,.  ♦  ♦ ;  °  “°te  slze  0  5 

higher  than  20  keV,  but  the  1  20..  —simulated 

transmission  under  20  keV  is  almost  H  10--  transmission, 

flat.  The  low  transmission  in  the  high  o - 1 - 1 - 1 - 1 - 1  1  mm’ no  waviness 

energy  range  can  be  explained  by  0  20  40  60  80  100 

the  non-linearity  of  the  lens.  The  Photon  energy  (keV) 

Outermost  channels  are  more  bent  Figure  49.  Transmission  of  the  central  part  of  the 


0  20  40  60 

Photon  energy  (keV) 

Figure  48.  Transmission  of  whole  tapered  lens. 


Pinhole  size:  1 


Pinhole  size:  0.5 


Simulated 

transmission, 

1  mm,  no  waviness 


Photon  energy  (keV) 

Figure  49.  Transmission  of  the  central  part  of  the 
tapered  optic. 
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(8.4  keV) 


up-down 


*960-1200 
H  720-960 
□  480-720 
*240-480 
*0-240 


right-left 


(35  keV) 

(fc-T  l\C  *  I 


Figure  51.  Two  dimensional  scans  of  the  output  of  the  tapered  optic. 

than  the  central  channels,  so  they  have  low  or  zero  transmission  at  the  higher  photon 
energies.  This  is  confirmed  by  the  higher  transmission  measurement  at  the  central  part 
of  the  lens,  shown  in  Figure  49  and  the  transmission  uniformity  measurement  at 
different  energies,  shown  in  Figure  51.  Transmission  uniformity  of  the  lens  was 
measured  by  putting  a  200  pm  pinhole  in  front  of  the  lens  and  performing  a  two 
dimensional  scan.  At  8  keV,  the  whole  lens  transmits,  although  we  can  see  the 
transmission  is  not  completely  flat  due  to  the  non-linearity  and  the  defects  of  the  lens, 
such  as  blocked  channels.  At  25  keV,  transmission  of  the  whole  lens  drops  to  25%. 
Most  of  the  lens  is  still  transmitting  at  25  keV,  but  the  outer  channels  have  lower 
transmission  than  the  central  part.  At  35  keV,  the  lens  has  a  narrower  transmitting 
region.  The  measurements  show  that  the  lens  has  undesirable  nonlinearity. 

Studies  were  performed  on  another  optic  to  investigate  the  distributions  of 
defects  which  diminish  transmission.  A  poorly  transmitting  optic  was  cut  into  two 
pieces.  The  transmission  of  each  piece  is  shown  in  figure  50.  The  transmission  of  the 

pieces  compared  to  the  original  is 
shown  in  figure  52.  This  work  is 
ongoing. 


Simulations  plotted  in  Figure 
47  and  Figure  48  as  the  solid  lines 
assume  when  the  lens  is  uniformly 
bent.  The  simulations  used  the 
same  roughness  and  waviness 
correction  as  that  used  for  fiber  D  in 
Figure  5.  A  linear  tapered  lens  is 
ideal  for  mammography,  because, 
as  shown  above,  any  curvature  can 
cause  transmission  reductions. 
However,  the  lens  currently 
available  has  a  nonlinear  profile. 

The  profile  is  not  circular,  but  to 
simplify  the  simulation,  we  assume 
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Figure  50 .  Transmission  of  two  pieces  of  a  tapered 
optic. 
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that  the  profile  is  circular.  A 
uniform  bending  will  provide  the 
best  transmission  performance  for 
a  given  bending  angle,  so  the 
simulation  will  overestimate  the 
lens  performance.  A  further 
simplification  assumes  the  lens 
has  a  fixed  channel  size,  although 
the  channels  actually  taper  with 
the  lens.  A  more  sophisticated 
lens  simulation  program  is  under 
development. 

The  geometry  for  uniformly 
bent  nonlinear  lenses  is  shown  in 
Figure  53.  The  simulation  method 
Figure  52.  Transmission  of  original  optic  compared  to  two  for  a  uniformly  bent  lens  is  an 

pieces.  extension  of  the  simulation  for 

uniformly  bent  polycapillary  fibers. 

The  simulation  is  done  by  sampling  channels  in  the  input  cross  section  of  the  lens  with 
an  even  step  size.  In  the  simulation  for  a  bent  fiber  it  is  assumed  that  the  capillary 
bends  in  -y  direction,  and  z  direction  is  along  the  input  axis  of  the  fiber  In  order  to  reuse 
the  old  simulation  code,  the  coordinate  of  the  simulation  x-y-z  is  rotated  to  x’-y’-z’, 
where  z’  is  the  direction  of  channel  axis  at  the  input  end  and  -y’  is  the  channel  bending 
direction.  Assume  the  source  position  and  the  capillary  position  at  input  end  is  (xs\  ys’, 
zs’)  and  (0,  0,  Zj’)  respectively,  then  the  source  capillary  distance  is  z,’-  zs’.  The  origin  of 
the  coordinate  is  selected  to  be  the  focus  point  of  the  lens. 

Another  parameter  needed  by  the  simulation  is  the  bending  curvature  of  the 
selected  channel.  This  value  varies  with  the  position  of  the  channel.  The  calculation  of 
the  bending  curvature  is  demonstrated  in  Figure  54. 

First  it  is  assumed  that  the  channels  are  evenly  distributed  on  both  input  and 
output  cross  section.  Let  input  and  output  diameter  of  the  lens  be  r}  and  r0,  and  focal 
distance  o7  be  f,  and  the  lens  length  7c  be  L,  then 
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Figure  53.  Geometry  for  Lens  simulation. 
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A  B  =  -yfx^  +  _y.2  ’  CE  =  AB  *  rD  /  r, 

f  /(/  +  L)=AB  /CD'  d  =  CE  -  CD 

For  a  small  bending  angle, 

cc*L/ R  =  d/(L/ 2) 

where  R  is  the  bending  radius,  d  is  the  distance  between  point  D  and  E  as  shown  in 
Figure  54,  and  a  is  the  bending  angle.  Thus 


A  L  C 


Figure  54.  Curvature  estimation  for  lens  simulation. 


R  =  L2  /  2d  (i) 

The  simulated  results  in  both  cases  show  higher  transmission  than  the 
experimental  results.  This  indicates  that  the  lens  has  more  reduction  in  transmission 
due  to  the  bending.  This  is  expected,  because  the  lens  does  not  have  a  uniform 
bending. 

Manufacturing  Issues 

The  only  practical  method  to  produce  larger  linear  tapered  anti  scatter  optics  is  to 
couple  several  smaller  optics  together  to  form  one  large  optical  component.  The 
method  we  propose  to  attach  the  individual  optical  elements  together  is  quite  similar  to 
the  method  that  is  used  by  the  company  to  align  polycapillary  fibers.  These  multifiber 
devices  consist  of  thousands  of  individual  polycapillary  fibers  that  are  strung  through 
rigid  grids  to  produce  a  single  large  optic.  The  individual  fibers  must  be  aligned  to  within 
a  few  microns  in  position  and  to  less  than  a  milliradian  in  angular  orientation.  These 
requirements  are  more  stringent  than  those  required  in  this  project.  The  alignment 
technology  is  not  likely  to  be  a  problem.  However,  first  high  quality  monolithic  tapers 
must  be  reproducibly  manufactured. 

XOS  has  developed  the  capability  to  make  long  tapers  with  diameters  up  to  5 
mm  with  transmission  which  indicate  the  future  potential,  but  which  is  currently  less 
than  half  the  expected  value.  Work  is  necessary  to  improve  profile  control  (linearity), 
increase  transmission  efficiency,  and  increase  achievable  optic  diameter.  Single 
examples  have  been  manufactured  which  demonstrates  sufficiently  good  performance 
to  support  the  medical  imaging  research  and  demonstrate  the  feasibility  of  the 
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approach.  However,  considerable  work  is  necessary  to  achieve  optics  suitable  for 
clinical  use. 

There  has  been  substantial  progress  in  achieving  the  necessary  manufacturing 
control  to  repeatably  produce  high  quality  optics.  The  manufacturing  technique  was 
carefully  studied,  the  critical  problems  were  identified,  and  a  plan  was  prepared  to 
modify  the  manufacturing  process.  The  first  step  was  to  increase  the  controllability  and 
reproducibility  of  the  manufacturing  process.  The  development  of  monitoring 
automation  and  feedback  systems  was  achieved.  As  a  result  of  the  improvement  in 
process  control,  systematic  and  rapid  improvement  has  been  made  in  the  development 
and  manufacturing  of  new  optics. 

Initially,  glass  preforms  and  prototype  optics  were  produced  in  Russia.  They 
were  of  relatively  poor  quality,  and  no  two  optics  were  equivalent  in  performance,  but 


Table  3 


PROBLEMS 

ACTION  NECESSARY 

•  decrease  in  open  area 

•  better  temperature  profile  control 

•  better  tension  control 

•  evaluate  tension  -  temperature 
trade  off 

•  inclusions,  primarily  glass  particles 

•  improved  handling  and  bundle 
assembly  procedures 

•  static  electricity  abatement 
including  controlling  for  humidity 

•  laminar  flow  hoods 

•  cast  (curvature) 

•  improved  alignment  of  feed 
mechanism,  furnace,  and  puller 

•  radially  symmetric  uniformity  of  heat 
zone 

•  longer  cooling  zone 

•  better  control  of  air  flow 

•  OD  control 

•  improved  temperature  and  tension 
control 

•  real  time  monitoring  of  OD  with 
closed  loop  feedback 

•  waviness  in  glass  surface 

•  elimination  of  internal  vibration 
sources,  e.g.,  friction,  motors  being 
operated  slower  than  designed, 
turbulent  forced  air  cooling 

•  isolate  from  external  vibration 
sources 

•  change  to  cutter  approach  which 
doesn’t  cause  large  periodic  shocks 

•  more  homogenous  glass  (optical 
quality) 
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they  served  as  proof-of-principle  devices.  Because  of  the  technical  and  political 
difficulties,  Russian  optics  and  preforms  are  no  longer  available.  When  it  became 
apparent  that  having  the  optics  manufactured  in  Russia  was  not  viable  long  term,  X- 
Ray  Optical  Systems,  Inc.  (XOS)  commenced  developing  a  U.S.  manufacturing 
capability. 

XOS  chose  to  work  with  the  polycapillary  subcontractor  that  has  the  best 
controlled  process  for  making  polycapillaries  which  work  well  at  8  keV.  The  glass  is 
formed  by  a  continuous,  steady  state,  vertical  pulling  process,  with  as  much  as  1,000 
meters  in  a  single  production  run.  Several  manufacturing  runs  were  conducted  with 
evaluations  of  the  polycapillaries  done  between  each  run.  Adjustments  to  the  process 
were  made  for  each  new  run.  XOS  personnel  participated  in  some  of  the  runs.  The 
polycapillary  performance  and  the  manufacturing  process  were  carefully  analyzed  to 
determine  what  changes  were  needed  in  manufacturing  to  obtain  the  desired 
polycapillary  capabilities.  This  manufacturer  does  not  use  computer  control.  There  was 
no  automatic  feedback  in  the  manufacturing  system.  The  system  was  controlled  by 
having  an  operator  observe  that  a  condition  such  as  the  outer  diameter  of  the 
polycapillary  or  the  furnace  temperature  had  drifted  from  the  desired  point  and  then 
manually  adjusting  one  of  the  control  variables.  The  polycapillary  performance  at  high 
energies  was  not  acceptable.  Problems  which  were  identified,  and  the  actions 
determined  to  be  necessary  are  identified  in  the  Table  3.  XOS  determined  that  this 
supplier  was  not  able  and/or  willing  to  make  the  necessary  equipment,  personnel,  and 
operational  changes.  With  this  information,  XOS  modified  an  in-house  glass  pulling 
system  to  investigate  the  effect  of  making  some  of  the  changes  which  were  determined 
to  be  necessary. 

A  systematic  evaluation  the  major  manufacturing  parameters  was  not  possible 
initially  as  the  first  prototypes  made  were  pulled  on  a  setup  where  there  was  little 
control  over  the  manufacturing  parameters,  and  little  ability  to  record  the  parameters 
during  the  forming  process.  The  initial  setup  was  put  together  inexpensively  to  provide 
some  indication  of  the  direction  to  proceed.  It  enabled  a  definition  of  a  general  range 
for  some  of  the  manufacturing  parameters,  such  as  pulling  speed,  so  that  acquired 
information  could  be  used  to  design  a  more  sophisticated  pulling  setup.  One  source  of 
variation  in  the  prototype  monolithic  optics  came  from  the  uncontrolled  air  flow  through 
the  furnace. 

Modeling  of  the  optic  shaping  process  was  conducted.  Modeling  provides  some 
understanding  of  the  underlying  relationships  between  the  manufacturing  parameters 
and  the  resulting  optic.  Constant  feedback  of  the  experimental  results  to  increase  the 
accuracy  of  the  model  is  an  important  input  to  the  scientific,  quality,  and  production 
methods.  The  modeling  is  based  on  a  one-dimensional  finite-differences  algorithm. 

The  preform  (large  outside  diameter  polycapillary  bundle)  is  split  into  finite  elements. 
This  approach  allows  the  thermal  and  mechanical  properties  of  each  segment  to  be 
calculated.  At  each  time  step,  the  ends  of  the  preform  are  moved  according  to  a 
predetermined,  but  likely  variable,  speed.  The  formation  of  the  lens  occurs  in  the 
heated  area  and  the  glass  will  flow  according  to  viscosities.  Since  the  total  volume  of 
each  element  is  kept  the  same,  the  elements  that  are  stretched  the  most  end  up  with 
the  smallest  diameter.  The  radiative  heating  by  the  furnace  as  well  as  radiative  cooling 
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Length  (mm) 


Figure  55.  Variation  in  optic  profile ,  compared  to  design. 

from  the  preform  must  be  considered.  The  thermal  profile  for  the  preform  can  be 
calculated  for  points  along  the  center  axis  of  the  furnace. 

A  second  generation  pulling  tower  was  constructed  and  instrumented  to  allow 
flexibility  in  controlling  different  parameters  and  to  record  the  actual  status  of 
parameters  during  the  pulling.  The  tower  uses  servo  motors  with  encoders  to  displace 
the  preform  and  furnace;  multiple  thermocouples  to  monitor  the  temperatures;  and  load 
cells  to  monitor  the  pulling  forces.  So  far  a  few  long  tapered  monolithic  optics  have 
been  formed  using  this  system.  There  are  numerous  independent  and  dependent 
variables  that  effect  the  quality  and  reproducibility  of  a  lens  design.  In  addition  to  the 
parameters  already  known  to  effect  the  reproducibility  and  quality  of  the  lens,  it  is 
expected  that  others  will  be  revealed.  An  important  part  of  the  project  will  be  to 
establish  the  relative  importance  of  each  parameter  through  a  systematic,  multivariable 
experimental  approach.  The  goal  being  to  seek  a  rapid  high-quality  path  to  large 
tapered  monolithic  polycapillary  lens  production.  Several  of  the  development  steps  will 
proceed  concurrently.  An  example  of  this  is  the  storage,  treatment,  and  holding 
methods  of  the  preforms  can  be  investigated  while  the  lens  formation  process  is 
proceeding  through  its  initial  data  collection  and  analysis  steps.  Although  the  ultimate 
quality  of  the  lens  is  dependent  on  each  of  these  paths,  the  synergism  is  likely  not 
realized  until  the  relative  significance  of  the  basics  is  established.  The  preform  material 
will  be  characterized  for  open  area,  x-ray  transmission,  uniformity  of  outside  diameter. 
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Figure  56.  Profiles  of  two  recently  manufactured  matched  linear  tapers. 

General  observations  of  the  presence  of  inclusions,  cleanliness,  and  uniformity  of  color 
will  be  noted.  Then  the  appropriate  method  of  storage  prior  to  use  needs  to  be 
investigated.  Precautions  to  protect  the  material  from  dust  and  moisture  will  be  taken. 

There  is  a  strong  interdependence  of  tension  on  the  preform,  temperature  in  the 
various  heating  zones,  speed  of  the  pull,  and  mechanical  feedback.  Extensive 
correlation  of  the  data  from  the  forming  process  with  the  lens  performance  will  be 
needed  to  optimize  the  system.  Additional  considerations  include  the  heat  capacity  of 
the  glass  and  time-temperature  response  of  the  furnace.  Uniformity  of  gas  flow  through 
and  around  the  glass  during  the  pulling  step  is  important  to  monitor.  To  assure 
uniformity  and  appropriate  use  of  heat  zones,  other  parameters  to  consider  include: 
baffling,  furnace  diameter,  and  correlation  of  furnace  motion  with  the  movement  (i.e. 
flow)  of  the  preform.  In  a  motor  driven  system,  there  are  vibrations,  motor  inertia,  slew 
rates,  as  well  as  bearing,  slide,  and  belt  qualities  that  must  be  evaluated.  Other 
mechanical  conditions  to  evaluate  include  how  the  preforms  are  held  and  supported  to 
maintain  uniform  alignment  and  minimize  random  and  systematic  perturbations.  Despite 
the  number  of  parameters  to  be  considered,  only  a  few  are  thought  to  be  critical  to  first 
order,  others  affect  the  quality  in  diminishing  amounts.  In  fact,  most  will  not  become  an 
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■  issue  of  concern  until  the  grosser  variable 

are  better  controlled.  Fairly  early  in  the 
process,  XOS  has  made  significant  progress 
in  the  design  and  construction  of  monolithic 
polycapillary  optics. 

The  first  example  is  a  set  of  17  keV 
focusing  optics  designed  for 
microfluorescence  analysis.  The  profile  of 
the  optic  was  measured  with  a  microscope- 
based  MicroVu  video  coordinate  measuring 
system  at  XOS.  The  precision  of  the  system 
is  about  5  pm.  Figure  55  shows  the 
measured  profile  compared  with  the  design 
profile.  Deviations  along  the  optic  axis 
between  the  two  profiles  are  also  shown  in 

Figure  57.  X-ray  transmission  image  of  Fi9ure  55'  More  recent|y'  **°  verV  we" 

linear  taper  shown  in  figure  56.  matched  magnifying  linear  tapers  were 

produced  for  placing  in  a  multi-taper  jig.  The 
profiles  of  the  two  optics  are  shown  in  Figure  56.  Although  the  optic  profiles  are 
excellent,  x-ray  transmission  efficiencies  of  all  four  optics  fabricated  so  far  are  pretty 
bad  (less  than  10%).  Based  on  x-ray  characterization  on  both  original  bundles  and  long 
linear  tapers,  we  conclude  that  the  poor  x-ray  transmission  is  due  to  randomly 
distributed  defects  in  the  original  bundles.  Figure  57  shows  an  x-ray  image  of  the  linear 
taper,  showing  some  areas  of  poor  transmission. 


Imaging  Measurements 
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Figure  58.  Contrast  enhancement  measured  for  a  small  scanned  capillary  optics  using  a  phantom 
consisting  of  a  Lucite  block  with  drilled  holes. 
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Figure  59.  Image  of  contrast  figure  60.  Image  of  contrast  11  lv'  ' ,v  ,VJ 

phantom  made  with  scanned  phantom  produced  with  conventional  grid  cases 

0pM-  conventional  grid.  are  shown  in  Fi9Ure  58. 

An  image  made  with  a  scanned  optic  and  with  a 
conventional  grid  are  shown  in  Figure  59  and  Figure  60.  Improvements  have  been 
made  in  the  scanning  system  to  reduce  image  artifacts.  Scan  lines  which  were  visible 
in  previous  images  are  not  apparent  in  Figure  59. 

The  optic  was  also  used  to  measure  the  resolution  enhancement  due  to 
magnification.  A  plot  of  the  measured  MTF  is  shown  in  Figure  61.  At  the  0.1  MTF 
level,  the  resolution  improved  from  4.5  Ip/mm2  without  the  optic  to  8.9  Ip/mm2  with  the 
optic.  The  MTF  of  the  very  linear  optic  from  Figure  56  is  shown  in  Figure  64.  It  shows 
nearly  ideal  performance  despite  the  low  transmission. 


Contrast 
enhancement 
measurements  have  been 
performed  with  a  small 
scale  (2  mm  diameter) 
optic  as  part  of  an  on¬ 
going  collaboration  with 
the  Univ.  of  Wisconsin.19 
The  geometry  for  the  scan 
is  shown  in  Figure  46,  and 
the  resulting  contrast 
enhancement  compared  to 
the  no  optic  and 
conventional  grid  cases 
are  shown  in  Figure  58. 


Artifact  Reduction 

Modulation  Transfer  Function  (MTF)  is  the  most  fundamental  measurement  of 
spatial  resolution  used  in  radiology.  The  standard  technique  is  to  image  a  slit, 
determine  the  line  spread  function  (LSF),  and  compute  the  Fourier  transform.  The  MTF 
measured  for  the  mammography  system  with  capillary  optics,  compared  to  that  without 


Modulation  Transfer  Funotion 
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the  capillary  optics,  will  give  the  resolution 
improvement  and  also  can  be  used  to 
diagnose  artifacts.  For  some  lenses,  the 
transmission  is  not  always  uniform,  as  shown 
in  figure  62,  which  is  an  image  of  the  output 
from  a  lens. 

The  structure  of  the  capillary  bundles 
was  obvious.  In  this  case,  the  LSF  could  not 
be  obtained  by  imaging  a  slit  without 
background  elimination.  The  background 
deducted  slit  image  is  the  slit  image  divided 
by  the  image  with  capillary  only.  However, 


Figure  61.  Modulation  transfer  function. 


this  involves  image  registration.  The 
information  for  background  that  can  be  used 


for  image  registration  is  too  little  in  the  slit 
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Figure  63.  Magnified  image  of  the 
background  with  lens  structure. 

image.  An  alternative  method  is 
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Figure  64.  MTF  of  newest,  most  linear  taper  at  different 
source  to  optic  distances.  MTF  is  nearly  ideal. 

calculating  the  LSF  as  the  derivation  of  the  edge 
spread  function  (ESF).  The  edge  image  left  enough 
background  for  registration  so  that  it  could  be  a 


feasible  method  in  our  case. 


Figure  62.  Magnified  image  of  a  A  registration  algorithm  (FMI-SPOMF)  proposed 

lead  edge.  Real  image  size:  by  Chen,  et  al. 20  was  used  in  this  work.  This  is  a 

5(mm)  X 5 (mm).  method  to  match  a  two-dimensional  image  to  a 

translated,  rotated  and  scaled  image.  The  approach 
consists  two  steps:  the  calculation  of  a  Fourier-Mellin  Invariant  (FMI)  descriptor  for  each 
image  to  be  matched,  and  then  matching  of  the  FMI  descriptors.  FMI  descriptors  are 
translation  invariant.  The  matching  of  the  FMI  descriptors  is  to  find  out  the  rotation  and 
scaling,  and  achieved  using  a  symmetric  phase-only  matched  filtering  (SPOMF).  In  our 
case,  there  is  no  change  in  scaling,  but  rotation  is  probably  involved.  When  the  rotation 
was  found,  image  translation  is  found  by  SPOMF  method. 

Images  with  and  without  edge  are  shown  in  Figure  62.  and  Figure  63.  The  edge 
was  made  of  a  lead  plate.  The  size  of  each  image  is  fifty  by  fifty  in  pixels.  They  were 
taken  in  an  experimental  digital  mammography  system  with  a  computed  radiography 
digital  phosphor  plate.  The  digital  CR  plate  has  limited  resolution,  5  Ip/mm,  and  is 
usually  considered  to  be  not  good  enough  for  clinical  mammography.  But  the  effective 
resolution  of  the  system  could  be  improved  by  the  magnification  of  the  capillary  optics. 
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Figure  66.  Reprojection  of  a  two- 
dimensional  edge  image  into  a  finely 
sampled  ESF. 


Figure  65  Three-dimensional  edge  image  after  background 
subtraction. 


The  magnification  of  the  capillary  lens  used  in  the 
system  is  1.86,  so  that  the  effective  resolution  should  be  improved  by  a  factor  of  1.86. 
The  registration  algorithm  was  implemented  in  IDL.  With  the  FMI-SPOMF  registration 
algorithm,  no  rotation  was  found.  The  translation  was  then  found  by  SPOMF  algorithm. 
The  difference  between  the  usual  correlation  method  and  SPOMF  method  is  that  it  only 
uses  the  phase  information.  The  phase-only  correlation  function  has  sharper  peak  than 
normal  correlation  function.  Sub-pixel  resolution  could  be  achieved  theoretically. 
However  it  was  found  that  the  resolution  is  not  good  enough  in  our  case.  This  may  be 
caused  by  the  relatively  small  image  size.  The  background-subtracted  result  was  further 
optimized  by  manual  shifts  in  sub-pixel  range.  The  result  after  manual  optimization  was 
shown  in  Figure  65  .  The  periodic  background  is  totally  gone.  A  cubic  spline 
interpolation  method  was  used  for  the  best  result  while  shifting  the  background  image 
relative  to  the  edge  image. 

The  presampling  MTF  includes  unsharpness  of  the  detector  and  the  sampling 
aperture.  In  order  to  eliminate  the  aliasing,  a  finely  sampled  ESF  is  obtained  with  a 
slightly  angulated  edge  in  a  single  exposure  as  illustrated  in  Figure  66.21,22  The  angle  in 
our  edge  image  in  Figure  65  was  calculated  to  be  around  6  °.  The  resultant  ESF  is 
plotted  in  Figure  67  .  A  direct  calculation  of  MTF  is  to  take  the  derivation  of  the  ESF 
and  result  a  LSF,  then  take  Fourier  transform.  As  it  can  been  seen,  the  ESF  in  our  case 
is  not  smooth  enough  for  a  direct  calculation.  An  alternative  technique  is  to  use  an  ESF 
fitting  procedure.23  In  the  fitting  method,  the  ESF  is  represented  with  a  term  (1-exp)  and 
an  error  function  (erf)  as  in  equation  (1).  MFT  is  calculated  by  equation  (2),  where  a,  b, 
c  and  d  are  fitting  parameters  from  equation  (1).  The  fitting  curve  is  plotted  in  Figure  67 
as  the  solid  line. 


Direct  Digital  Mammography  Using  Capillary  Optics 


39 


ESF(x)  =  a{  1  -  exp(-Z>  |  x  -  x0 1)}  +  c  ■  erf  ( du 2  \  x  -  x0 1)  0) 

mF{f)  =  c-exp(-,r2/2/rf)  +  «(l  +  4;r2/2/62)-'  (2) 

(c  +  a) 


Figure  67  Edge  spread function  and  its  fitting 
curve. 
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Figure  68  Calculated  MTF  compared  with  the 
MTFs  with  scan  optics,  and  without  optics. 


The  calculated  MTF  was  compared  with 
that  from  scanned  optics  in  Figure  68  .  The 
MTF  from  scanned  optics  was  calculated 
from  a  slit  image.  Optic  structure  was 
smeared  out  with  scanning,  so  no 
background  deduction  is  necessary.  An 
ideal  MTF  was  also  plotted  in  Figure  68  . 

An  ideal  MTF  is  the  one  when  we  assume 
the  lens  does  not  introduce  only  the  1 .86 
factor  magnification  and  no  image 
degradation.  It  is  measured  with  a  slit 
image  without  optics,  and  its  result  was 
multiplied  by  the  magnification  of  the 
capillary  lens.  The  measured  MTF  with 
lens  should  no  better  than  the  ideal  MTF. 
However,  for  some  reason,  the  MTF  for 
stationary  lens  is  better  than  the  ideal  MTF. 
It  is  possible  that  the  background 
elimination  process  failed  the  MTF 
measurement.  Another  possibility  is  that 
the  magnification  of  the  lens  is  different 
from  1 .86,  since  the  magnification  was 
measured  separately  at  a  different  time. 
More  experiments,  necessary  to  determine 
the  cause  of  the  disagreement,  are 
underway. 

Demaqnifvinq.  Focusing  Monolithic  Optics 
Monolithic  optics  similar  to  the 
collimating  optic  shown  in  can  be  used  to 
minify  an  x-ray  beam.  Some 
characterization  of  such  optics  has  been 
performed.  The  first  measurement 
characterized  an  optic  using  a  point 
source,  a  micro-focused  transmission  x-ray 
source  (spot  diameter  about  5  |am) 
modified  from  an  old  SEM.  Anode 
materials  made  of  Cu  and  Mo  are  coated 
on  a  thin  Be  vacuum  window.  An  energy- 
sensitive  Amptek  detector  with  an  energy 
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Figure  70.  Focal  spot  scans  of  focusing  optic. 


Amptek  detector 


t 


X-Y-Z  stages 


Figure  69.  Schematic  of  the  experimental  setup  for  optic  characterization. 

resolution  of  250  eV  at  5.9  keV  was  used  for  x-ray  detection.  The  crystal  size  of  the 
detector  is  7  mm2.  The  detector  was  placed  roughly  at  the  output  focal  plane  of  the 
optic,  8.5  mm  away  from  the  exit  of  the  optic,  to  allow  the  detection  of  the  output  beam 
from  the  entire  optic.  The  polycapillary  focusing  optic  and  detector  were  mounted  on  x- 
y-z  stages  which  allows  precise  alignments  among  the  source,  optic,  and  detector. 
Figure  69  shows  the  schematic  of  the  experimental  setup.  The  alignment  of  the  optic 


Figure  71.  Knife  edge  scan  at  a  knife  edge  to 
source  scan  distance  of  103.5  (*)  and  123.5  mm 
(B). 


input  focus  to  the  source  is  achieved  by 
moving  the  optic  in  all  three  orthogonal 
directions  until  the  transmitted  flux  is 
maximized. 

The  output  focal  spot  of  the  prototype 
optic  was  at  Cu  Ka  and  Mo  Ka.  The 
measurement  was  performed  by  scanning  a 
tungsten  knife  edge  plate  across  the  focused 
beam  on  the  focal  plane.  The  first  derivative 
of  the  measured  intensity  as  a  function  of  the 
position  of  the  knife  edge  gives  the  intensity 
distribution  of  the  focused  beam.  This  is 
shown  in  Figure  70,  with  the  simulation 
results  as  comparisons.  Good  agreement 
between  the  measurements  and  the 
simulations  indicate  good  alignment  of  the 
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Figure  73.  Derivative  of  knife  edge  scan  for 
optic  to  knife  edge  distances  of  103.5  mm  (•) 


knife  edge  -  optic  distance  (rrm) 


and  123.5  mm  (R). 


Figure  72.  Width  of  derivative  of  knife-edge 
scan  versus  distance  to  optic  for  monolithic 
optic. 


capillaries  toward  the  output  focus  and 
that  the  simulation  can  reliably  predict 
the  divergence  of  the  output  beams  from 

the  capillaries.  The  measured  spot  sizes  (FWHM)  are  44  pm  for  Cu  Ka  and  21  pm  for 
Mo  Ka,  respectively. 


A  second  set  of  measurements  were  performed  at  beamline  X23A3  at  the 
National  Synchrotron  Light  Source  (NSLS)  at  Brookhaven  National  Laboratory  (BNL). 
This  is  a  white  beam  bending  magnet  line  with  a  maximum  white  beam  flux  of 
approximately  13  W/cm2.  The  size  of  the  beam  could  be  adjusted  by  two  sets  of 
collimating  slits.  The  optic  was  mounted  on  a  precision  motion  controlled  four  axis  stage 
(up-down,  left-right,  tilt,  and  rotate)..  A  knife  edge  or  pinhole,  mounted  on  another 
precision  stage  (left-right,  forward-backward),  could  be  placed  between  optic  and  the 
ion  chamber.  A  small  portable  laser  was  used  for  rough  alignment  of  the  optics  with  the 
beam. 

The  spot  size  and  also  the  focal  distance  were  measured  using  a  knife 
edge  that  was  scanned  horizontally  across  the  beam.  Two  scans,  one  near  to  the  focal 


Table  4.  Results  for  monolithic  focusing  optic. 


X-Ray 

Energy 

(keV) 

Spot  size 
(mm) 

Trans¬ 

mission 

(%) 

measured 
Gain  for 

350  pm 
pinhole 

calculated 

Gain 

350  pm 
pinhole 

calculated 

Gain 

90  pm 
pinhole 

calculated 

Gain 

10  pm 
pinhole 

6 

0.09 

36 

78 

81 

645 

911 

8 

0.08 

49 

96 

110 

933 

1359 

10 

0.09 

39 

83 

87 

624 

842 

12 

0.09 

39 

74 

87 

654 

903 

white 

0.17 
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11 

89 
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Figure  74.  Spectrum  from  Am  241  taken  with  CZT  single 
pixel  detector. 


spot  of  the  multifiber  optic  and  one 
20  mm  closer  to  the  optic,  are 
shown  in  Figure  71.  The 
derivatives  of  these  scans  are 
shown  in  Figure  73.  The  scan 
close  to  the  focal  spot  has  a 
smaller  width  than  the  scan  far 
away  from  the  focal  spot.  A  plot  of 
the  width  of  the  scan  derivative  as 
a  function  of  distance  from  the 
multifiber  optic  is  given  in  Figure 
72.  Resulting  calculated  and 
measured  gains  are  shown  in  table 
4 


Direct  Digital  Detectors 

CZT  Detector 

Digirad,  Inc.  (formerly  Aurora  Technologies)  was  somewhat  delayed  by  the  late 
start  date  of  their  subcontract  with  the  University  at  Albany  due  to  contract  formalities. 
However,  a  1  cm  long  CZT  array  with  50  urn  pixels  has  been  manufactured  and  tested. 
The  interconnect  options  were  researched  and  the  decision  was  reached  to  employ 
indium  bump  technology  to  bond  the  detector  to  a  sapphire  or  quartz  interconnect 
board,  which  is  then  bonded  to  the  readout  chips.  This  approach  eliminates  the  need  to 
wirebond  to  the  detector,  and  also  reduces  the  pixel  size  and  consequently  the  leakage 
current.  The  first  prototype  was  built,  which  demonstrated  some  mechanical  instability. 
It  is  believed  this  is  fixable  with  epoxy.  The  detector  is  currently  at  Digirad  under 
investigation. 

Measurements  of  room  temperature  energy  bandwidth,  detector  quantum 

efficiency,  and  response 


0  5  10  15  20  25  30  35  40 

keV 


uniformity  were  performed  with  a 
single  pixel  bulk  CZT  detector. 
Measurements  of  pixel  to  pixel 
variation  using  a  scanned  slit 
were  performed  with  a  linear  CZT 
array. 

Single  pixel  CZT  detector 

A  single  pixel  detector  5.0 
mm  x  5.0  mm  in  area  and  1 .5  mm 
in  thickness  was  used  for  the 
measurements.  The  external 
aluminum  window  was  removed 


Figure  75.  Spectrum  of  copper  x-ray  source,  taken  with  single  ®  S  ec  OT  quan  um 

pixel  CZT  detector.  The  lower  level  discriminator  was  set  at  ©fficiency  measurement. 
the  lower  end  of  the  region  of  interest  (ROI)  marked  in  figure  The  spectrum  produced  by 

5.  the  CZT  single  pixel  detector  with 
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Figure  77.  Spectrum  of  copper  x-ray  source,  mm 

taken  with  HpGe  detector.  Lines  at  10.5  and 

12. 6  keV  are  lead  fluorescence  from  shielding.  Figure  76.  Pinhole  scan  of  CZT  single  pixel  with  band 
The  region  of  interest  (ROI)  is  marked.  of  contact  material. 

a  Am-241  radioactive  source  is 

shown  in  Figure  74.  The  spectrum  was  obtained  with  a  Tennelec  TC  241  linear 
preamplifier  connected  to  a  PCA-Multiport  multichannel  analyzer.  The  noise  threshold 
at  room  temperature  with  this  preamplifier  is  10  keV.  The  FWHM  of  the  Am  59.54  keV 
gamma  is  9.5  keV  at  room  temperature.  The  17  and  24  keV  lines  of  Am  could  not  be 
resolved. 

The  detector  quantum  efficiency  was  measured  in  comparison  to  a  high  purity 
germanium  (HpGe)  detector  with  a  13  mm  thick,  36  mm  diameter  crystal.  For  the  HpGe 
detector,  the  efficiency  was  considered  to  be  equal  to  unity  in  the  energy  range  in  which 
the  experiment  was  performed.  This  detector  was  connected  to  the  same  linear 
amplifier  and  multichannel  analyzer  as  the  CZT  detector.  The  distance  between  the 
source  and  the  100  pm  pinhole  was  50  cm.  The  distance  was  chosen  to  maintain  dead 
times  less  than  10%.  The  distance 


Figure  78.  Spectrum  from  radioactive  iodine 
taken  with  HpGe  detector. 


Figure  79.  Spectrum  from  radioactive  iodine  taken 
with  CZT  detector 
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between  the  plane  of  the  pinhole  and  the  detectors  was  kept  equal  to  5  mm  for  both 
detectors,  which  was  as  small  as  possible  given  the  HpGe  detector  geometry.  The 
small  separation  reduces  the  effects  of  scattered  or  fluorescent  radiation  and  ensures 
that  the  entire  beam  exiting  the  pinhole  is  intercepted  by  the  detector.  The  HpGe 
detector  was  calibrated  in  energy  using  the  Am  gamma  and  the  5.9  keV  line  of  Fe  55. 
The  CZT  detector  was  calibrated  using  the  Am  59.54  keV  line  and  a  pulser  set  to  equal 
and  to  half-height  of  this  line. 

Two  different  x-ray  sources  were  used  for  this  measurement.  First  was  a  copper 
x-ray  tube  with  tube  voltage  of  30  kV,  with  the  spectrum,  obtained  with  the  HpGe 
detector,  shown  in  Figure  77.  The  resultant  CZT  spectrum  for  the  same  count  time  is 
shown  in  Figure  75.  Interpretation  of  the  results  is  complicated  by  two  factors.  First, 
the  noise  threshold,  and  therefore  the  lower  level  discriminator,  is  above  the  peak  of  the 
spectrum,  so  only  the  higher  energy  bremsstrahlung  is  recorded.  Second, 
monoenergenic  high  energy  photons,  due  to  hole  tailing  effects,  contribute  to  counts  in 
a  long  Laurencian  tail  to  lower  energies,  as  seen  for  the  Am  gamma  in  Figure  74. 
Therefore  counts  due  to  photons  with  energy  above  the  threshold  are  lost  into  lower 
energy  channels.  The  detected  quantum  efficiency  for  this  source,  determined  roughly 
as  the  ratio  of  integrated  counts  in  the  regions  of  interest  shown  in  Figure  77  and  Figure 
75,  was  found  to  be  44  %.  A  similar  measurement  was  performed  for  a  radioactive 
iodine  source,  with  the  spectra  shown  in  Figure  78  and  Figure  79.  The  DQE  for  this 
source  is  42%.  An  imaging  detector,  which  counts  events  independent  of  energy, 
would  be  affected  by  the  hole-tailing  only  to  the  extent  that  the  ends  of  the  tail  are 
below  the  noise  threshold. 


mm 


mm 


Figure  80.  Pinhole  scan  of  CZT  single  pixel. 


Area  uniformity  measurements 
were  performed  on  a  bulk  single  pixel 
detector  by  scanning  the  active  area  of 
the  detector  with  a  300  pm  pinhole.  The 
result  for  the  first  detector  tested  is 
shown  in  Figure  76.  The  diagonal  band 
is  a  strip  of  contact  material  normally 
unimportant  for  the  operation  of  this 
detector  in  the  infrared.  The  scan  of  a 
second  detector,  produced  without  this 
contact  material,  is  shown  in  Figure  80. 
The  average  count  rate  in  the  3.5  mm 
central  area  was  normalized  to  unity. 
The  variance  in  that  region  is  8%,  3%  of 
which  is  due  to  Poisson  statistics. 

Some  part  of  the  variance  (4%)  is  due 


to  the  external  aluminum  window.  The  semiconductor  sensitivity  is  quite  uniform  over 
the  entire  pixel. 

Linear  array 
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Pixel  number  Pixel  number 

Figure  81.  Image  of  slit  taken  with  linear  CZT  Figure  82.  Background  associated  with  image  in 

array,  without  background  subtraction.  figure  1 1. 

The  CZT  linear  array  was 

attached  to  the  read-out  electronics  by  means  of  indium  bump  bonding  technology. 
Measurements  were  performed  only  for  the  central  pixels  because  mechanical  stresses 
caused  bonding  failure  for  the  end  pixels.  This  will  be  corrected  for  future  arrays  with 
an  addition  of  epoxy  for  mechanical  adhesion.  To  check  the  uniformity  of  the  response 
of  the  pixels  to  x-ray  radiation,  a  0.5  mm  slit  was  placed  between  the  detector  and  the 
copper  x-ray  source.  The  array  was  then  mounted  on  a  translation  stage  that  allowed 
movement  parallel  to  the  plane  of  the  slit.  The  result  without  background  subtraction  is 
shown  in  Figure  81 .  The  background,  counted  for  the  same  period  without  the  source 
on,  is  shown  in  Figure  82.  The  background  subtracted  result  is  shown  in  Figure  84. 
Moving  the  detector  to  the  left  causes  the  peak  to  move  to  the  right,  as  shown  in  Figure 
83.  Comparing  the  peak  heights  of  six  positions  along  the  central  part  of  the  array 
gives  a  peak  height  variance  of  25%.  To  establish  that  the  variance  was  due  to 
differences  in  pixel  sensitivity,  a  full  field  measurement  was  taken  with  no  slit.  Dividing 


Pixel  number  pixe|  number 


Figure  84.  Slit  image  from  Figure  86  with 
background from  Figure  85  subtracted. 


Figure  83.  Background  subtracted  image  of  slit  after 
moving  detector  to  the  left. 
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Figure  85.  CID  image  of  three  grids.  The  largest  has  open 
spaces  560  microns  wide  and  300  micron  wires.  The  smallest, 
unresolved  at  upper  right,  has  25  micron  holes.  The  grid  at 
center  has  50  micron  holes.  The  image  was  taken  with 
magnification  1. 02  onto  the  chip,  which  has  28  pm  pixels. 


the  background  subtracted  slit 
image  by  the  background 
subtracted  full  field  image 
causes  the  two  peak  heights 
shown  in  Figure  84  and  Figure 
83  to  agree  within  6%. 

Unfortunately,  problems 
have  developed  in  advancing 
this  technology  from  linear  to 
two  dimensional  arrays. 
Although  only  one-dimensional 
arrays  were  proposed  in  this 
project,  it  was  recognized  that 
the  eventual  outcome  should 
be  a  2D  array.  The  current 
technology  for  2D  CZT  arrays 
do  not  lend  themselves  to  the 
small  pixel  sizes  and  relatively 
low  x-ray  energies  (compared 
to  nuclear  gamma  imaging) 


required  for  mammography. 

CID  Detector 

Because  of  the  questions  in  development  of  the  CZT  technology  we  have 
investigated  a  second  technology,  Charge  Injection  Devices  (CID),  as  well.  Like  the 
more  common  CCD  (Charge  couple  device)  CID  technology  allows  for  rapid  array 
imaging.  Unlike  the  CCD,  the  CID  is  very  radiation  resistant.  CID  pixels  can  be  read 
repeatedly  and  nondestructively,  leading  to  very  high  dynamic  range.  CID  well  depths 
are  also  deeper  than  conventional  CCDs,  leading  to  a  substantial  reduction  in  blurring 
from  high  energy  x-ray  photons.  The  CID  chips  tested,  which  had  250x250  28  pm 

pixels,  were  too  thin  to  have 
good  x-ray  quantum  efficiency 
without  phosphor.  However, 
large  area  chip  technology  is 
already  available.  The 
measured  DQE  for  the  device 
was  23%  at  20  keV.  At  that 
enery  only  29%  of  the  photons 
are  expected  to  be  absorbed  in 
the  8  pm  thick  active  region. 
Detector  designs  with  thicker 
active  areas  will  have  larger 
DQEs.  We  are  working  with 
the  manufacturer,  CIDTec,  to 

Figure  86.  CID  image  ofLucite phantom  shown  in  Figure  42.  Study  the  application  of  their 
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new  x-ray  dental  imagers  (with  and  without  phosphor  coatings)  to  mammography. 
Radiation  damage  testing  measurements  have  been  completed  on  prototype  imagers. 
Initial  testing  of  a  prototype  imager  with  phospor  resulted  in  mechanical  stress  damage 
which  is  being  investigated.  A  second  prototype  for  direct  detection,  without  phosphor, 
will  be  delivered  in  the  next  few  weeks. 

Conclusions 

This  is  a  high  risk  project  with  potentially  very  high  gain.  Any  such  project  will 
have  inevitable  delays,  frustrations  and  changes  of  direction.  However,  the  delays  and 
problems  associated  with  this  project  have  been  fairly  minor.  Some  manufacturing 
problems  have  occurred  with  both  optics  and  detectors.  These  problems  are  being 
addressed  by  our  studies  of  defect  properties,  and  by  changes  being  adopted  in  the 
manufacturing  processes.  Notwithstanding  the  problems,  the  optics  tested  have 
performed  extremely  well.  A  number  of  promising  geometries  are  being  studied: 
collimating  optics  with  long  and  short  focal  lengths,  with  and  without  antiscatter  optics, 
monolithic  linear  magnifying  tapers,  and  monolithic  focusing,  demagnifying  optics.  The 
collimating  optics  have  transmissions  in  excess  of  30%  at  their  design  energies,  with 
collection  angles  ranging  from  8  to  12  degrees.  Scatter  rejection  is  very  high  from  all 
the  optics.  The  linear  tapers  resulted  in  nearly  ideal  contrast  enhancements  and 
simultaneously  increase  in  MTF  at  all  spatial  frequencies.  This  will  result  in  an 
improvement  in  resolution  regardless  the  inherent  spatial  resolution  of  the  detector. 

Just  as  important,  if  less  spectacular,  is  the  rapidly  growing  modeling  capability.  The 
modeling  has  lead  to  a  real  development  of  understanding  of  the  nature  of  polycapillary 
defects,  which  is  already  leading  to  improvements  in  the  manufacturing  processes. 
Further,  the  increased  confidence  in  the  modeling  allows  future  lens  geometries  to  be 
effectively  “tested”  in  simulation  so  that  design  parameters  can  be  rapidly  optimized. 

Direct  digital  detector  development  is  also  promising.  Two  diverse  detector 
technologies  have  been  investigated.  Both  show  considerable  promise  for  providing 
simultaneous  high  resolution  and  high  quantum  efficiency.  One  dimensional  “imaging” 
was  demonstrated  with  the  CZT  linear  array.  The  CZT  detector  has  higher  quantum 
efficiency  without  phosphor,  but  is  more  difficult  to  produce  as  a  two-dimensional  array 
with  sufficiently  low  read-out  noise  and  small  enough  pixels.  However,  this  is  largely  a 
problem  with  the  read-out  electronics.  As  better  electronics  are  becoming  rapidly 
available,  this  technology  becomes  more  attractive.  The  CID  technology,  while  it  would 
require  some  modification  for  phosphor  free  operation,  is  moving  rapidly  into  the  x-ray 
market,  particularly  in  dental  imaging.  CID  two  dimensional  array  technology  with  small 
pixel  sizes  already  exists.  This  is  another  promising  new  technology  for 
mammography. 
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